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ABSTRACT 

Velocity and temperature fluctuations from multi-level 
measurements in the marine boundary layer were examined to 
determine the validity of present formulations and prediction 
techniques. Spectral distributions were analyzed using an 
analog analyzer as opposed to the digital computer methods 
previously utilized at the Naval Postgraduate School. 

In general, the height variations of the dissipation 
rate of turbulent kinetic energy, momentum flux, and the 
temperature structure parameter observed N in this study agree 
with the overland expressions developed by Wyngaard. Most 
deviations from expected results were attributed to apparent 
anomalies in spectral results. 

It has been shown that the analog spectral analysis 
performed in this study produced information which is as 
reliable as that from digital methods. 
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I , INTRODUCTION 



This is a study of velocity and temperature fluctuations 
in the marine boundary layer. The data examined were obtained 
from multi-level measurements with tungsten wire velocity 
probes and platinum wire temperature probes. The measure- 
ments were made from the Naval Postgraduate School research 
vessel, ACANIA, Mean wind and temperature measurements are 
also utilized to complement and analyze the fluctuation data, 

A separate investigation on the mean data has been completed 
by Smedley (1975). 

The purpose of this study is to examine properties of 
small scale temperature and velocity fluctuations in the 
marine environment, on the basis of spectral distributions. 

A significant aspect of the study is that spectral distribu- 
tions were obtained utilizing an analog spectrum analyzer, 
rather than digital computer methods. The latter has been 
the method used previously in atmospheric turbulence experi- 
ments at the Naval Postgraduate School. The study includes 
descriptions of the variation with height of (1) e (the viscous 
dissipation of turbulent kinetic energy) and (2) (the 
temperature structure parameter). 
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II, THEORETICAL CONSIDERATIONS 



A . BACKGROUND 

This study is one of a continuing investigation of 
multiple-level overwater observations from a ship. The ship 
is the R/V ACANIA, shown in Figure 1, operated by the Depart- 
ment of Oceanography, Naval Postgraduate School. Because 
these were shipboard measurements, only high frequency small 
scale fluctuations could be measured. The following paragraphs 
will review the theoretical considerations used to interpret 
the small fluctuation data, Parallel reviews of this theory 
have been provided by Johnston (1974) and Bone (1974). The 
latter involved initial examinations of data obtained with 
the same instrumentation used in this study. 



B. TURBULENT KINETIC ENERGY BALANCE EXPRESSION 

In a simplified form, the turbulent kinetic energy 
balance expression for neutral conditions can be written 



3u 

3z 



= e 



( 1 ) 



where u represents the mean wind at height z and the term 
is the wind shear, Equation (1) states that mechanical 

r\' 

production of turbulent kinetic energy, -u ' w ' equals the 

molecular dissipation of turbulent kenetic energy, e. 



Equation (1) is often written as 

u’J*. e 

* 3z 



( 2 ) 
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Figure 1. The R/V AC AM A 
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where u* = -u'w' , and -u'w' represents the vertical transfer 
of horizontal momentum by turbulent fluctuations. 

For neutral conditions in a constant stress boundary layer, 
the following expression is valid for the wind gradient 



9u U * 



(3) 



3z Kz 

where K is Von Karman's constant. Combining Equations (2) 
and (3) yields an expression to estimate u on the basis of 
measurements of e at a known height z. 



= eKz 



Equation Q4) can be rewritten 

u * 

Ln e = -Ln z + Ln 



(4) 



(5) 



which is a form used in later analysis of the variation of 
£ with height . 

For non-neutral conditions the turbulent kinetic energy 
balance expression (Garratt 1972) is 

( 6 ) 



u * = eKz M — ) 



For an unstable case (R. < 0), 

v 1 ' > 



<f> (S) = (1 - 16 C) _l 

m 



(7) 



and £ = R. . 



For the stable case (R^ > 0), 
<t> m U) = 1 + 5 £ 



and 



S = 



R. 

l 

1 - 5 R. 



( 8 ) 
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where R. is the Richardson number. 

1 ' 



R i " 8/T OT/ 3z ) / (3 u/ 3z ) 2 



Equation (5) can be rewritten as 



Ln e = -Ln ( 2 / 



( 4> m (0 -£) 



u * 

) + Ln r— 
k 



(9) 



A plot of Ln e versus the first term on the right hand 
side of Equations (5) or (9) would yield a straight line with 
-1 slope. Figure 2 illustrates the possible effects of 
instability on the predicted plot of Ln z versus Ln z, where 
Ln <j> 2 (R^) is an empirical function of the Richardson number, 
Deviations from this relationship could be interpreted as 
either due to deviations from the constant flux assumption 
arising from the influence of wind-wave coupling or non- 
neutral stratifications. 

From the above discussion, e can be used to estimate u + 
values.. u + is an important parameter in boundary layer studies 
since it defines the intensity of turbulent transfer of 
momentum toward the surface. It is assumed to be constant 
with height in the first 30 meters. 

C. SPECTRAL SIMILARITY THEORY 

Kolmogorov's second hypothesis is the basis for estimating 
e and C^, (the temperature structure function parameter) 
from spectra of velocity and temperature fluctuations. It 
is based on similarity expressions for the inertial subrange, 
which is the wave number region where energy is transferred 
v/ithout dissipation from lower to higher wave numbers. 
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Ln Z 



Figure 2, An Illustration of the Effect of Stability on 
Variation of e with Height, 
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The following relationships were proposed by Kolmogorov 



for spectral densities in the inertial subrange, 




( 10 ) 



, - 5/3 

4> T (k) = c 2 c* k ' 



( 11 ) 



where <f>(k) and <J>,j,(k) are the one dimenensional spectra for 
velocity and temperature respectively and k is the streamwise 
component of wave number as defined below. 

e has already been discussed in it's relationship to u^, 
Another turbulence parameter of interest in this study is the 
temperature-structure parameter, c£, which appears in 
equation (11). It is important because it is related to the 
refractive-index-structure parameter, C^, by the expression 



is a primary parameter in the characterization of optical 
propagation in the environment. 

Height variations of are dependent on stability, as 
was the case for e. Predicted height variations for various 
stability conditions are shown in Figure 3. Wyngaard, et . al . 

^ 2/3 

(1971), noted that decreased with height as Z” ' for 

~ 4 / 3 

stable conditions and as Z ' for unstable conditions. 

Since velocity and temperature fluctuations are measured 
at a fixed point in the flow, the resultant spectra are 
realized at a "temporal" frequency, f, To obtain an e and 
C,p, the time (f) and. space (k) scales must be related by 
Taylor's (1938) "frozen turbulence" hypothesis, i.e., 




( 12 ) 
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. Ln Z 



Figure 3 




An illustration of the effect of stability on 
the variation of Ci, with height. 
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k 



2nf 



u 

where u is the mean wind speed at the measurement level. The 
term "frozen" implies the turbulence pattern remains unchanged 
as it sweeps past the probe. 

Equations (11) and (12) can also be written as 

f<J>(f) = k<Kk) = c e 2 / 3 k~ 2 / 3 (13) 

f$ T (f) = ^4> t ( k) = c 2 c 2 k~ 2/3 (14) 

where 4>(f) and <f>,p(f) are spectral density values with units 
of (m/sec) 2 /Hz and °C 2 /Hz respectively, 

Final forms of Equations (13) and (14), in terms of <j>(f) 
and <j)^,(f), which are measured, are 

/ v 2 / 3 \ 3 / 2 

e = (15) 

C 2 = k m f <j> ( f ) (16) 

c 

2 

where, empirically, c x = 0.5 and c 2 = 0.25. 

From Equations (15) and (16), with measured f, 4>(f)» 4>-p(f), 
and u, e and C^, can be determined at each level of interest. 
These were the final forms used to determine the variation 
with height of e and C 2 , values. 
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Ill, DATA COLLECTION AND SENSOR CALIBRATION 



A . THE PLATFORM 

Observations were obtained aboard the R/V ACANIA anchored 
off Monterey, California in Monterey Bay, Measurements were 
made at four levels on two masts spatially mounted on the 
forward deck of the ship. Figure 4 illustrates the sensor 
locations on these masts. Figure 5 shows the vane and probe 
arrangement, The vane maintained a one dimensional profile 
of the wind during measurements, 

B. INSTRUMENTATION 

The velocity fluctuation sensors were TSI model 1210 
probes with sliding support shields and tungsten wires. The 
shields permitted isolation of the sensing region for deter- 
mination of the undisturbed velocity, V , before and after 
each experiment. The platinum coated tungsten wire was small 
enough to resolve the viscous dissipation scale without any 
need to apply wire length corrections. The wire was fitted 
with plated ends for isolating the sensing area and thereby 
minimizing flow disturbance. 

Electronics associated with each probe were a TSI Model 
1054B Linearized Anemometer , a TSI Model 1056 Variable Decade 
Module and a Model 1057 Signal Conditioner. The signal 
conditioner was only used during the March 1974 observations. 
The anemometer had a linear frequency response from DC to 10 



19 



legend 



U Cup Anemometer 
T Quartz Thermometer 
q Humid iometer 
u 1 Hot Wire 
T 1 Platinum Wire 



Level 4 



Level 3 
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Figure 4. Mounting arrangements, 
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Figure 5 . Vane and probe arrangement . 



KHz and the variable decade module operated with a 0-60 ohm 
range. The signal conditioner had a linear frequency response 
from DC to 400 KHz and permitted voltage suppression in one 
volt steps from 0-29 volts. 

Sensor placement required exceptionally long cables, but 
any decrease in system frequency response had little effect 
in the frequency band of interest. 

In general the performance of the hot-wire system was 
excellent throughout the experiments. 

Temperature fluctuations were also measured with a TSI 
Model 1210 configuration using platinum wire sensors. 
Electronics used were the GTE Sylvania Model 140 thermosonde 
temperature modules, which can be used to make both single 
point temperature fluctuation and AT measurements. The 
temperature module is designed for measurements of temperature 
variations down to 0,001°C at frequencies up to 1 KHz. 

Mean measurement systems are described in Smedley's 
(1975) study. 

C. SENSOR CALIBRATION 

1 . Hot Wires Calibration 

In-situ calibrations of the velocity sensors were 
accomplished several times during each experiment. In this 
procedure, recordings v/ere made of both the cup anemometer 
RPM count and the corresponding hot wire voltage output. The 
sensor wind speed and voltage relationship is given by 

v 2 = a(u) 2 + b (17) 
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where v is the voltage and u is mean wind speed for a partic- 
ular level. The constants a and b are calibration curve slope 
and intercept respectively obtained from the in-situ calibra- 
tion. Figure 6 is a sample calibration curve constructed 
from in-situ values in one period. 

To convert the power spectral density levels to 
velocity units required a calibration factor such that 



c v* = u' (18) 

cm / s © c 

where c is the calibration factor in — L -r r ■ — , v' is the 

volt ’ 

voltage flucuation and u' is the velocity fluctuation. 
Differentiating Equation (17), we get 



( 



4v(u)* 

a 



) 



V 1 = u f 



(19) 



or 

c = (2 

To obtain the calibration factor, c, Equation (17) 
was initially solved for v using the appropriate sensor 
level's mean wind speed u and also the particular sensor's 
calibration curve slope and intercept. 

2 . Platinum Wires Calibration 

Wire resistance is related to temperature by the 



equation 

R = Ro(l + a(T - To)) (21) 

where R = resistance, T = calibration temperature, To = 0°C, 
Ro = resistance at To, and a = ,0039/°C (temperature 
coefficient of resistance for platinum). Each wire's 
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Figure 6. Hot wire calibration plot. 
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resistance was measured in the laboratory at a temperature T, 
to determine Ro, 



We can relate temperature fluctuations to wire 



resistance fluctuations by differentiating Equation (21), 



Recorded voltage fluctuations can then be related to 



resistance changes by the transfer function, H, in the 
expression 



where H was determined for the four bridges by varying the 
resistance of a TSI 1056 Variable Decade Probe Resistance 
and recording the corresponding voltage changes. Figure 7 
is a calibration curve of voltage versus resistance, where 
the slope determines the transfer function, H. 



temperature fluctuations substitute Equation (23) into 
Equation (22), giving 



dR = a Ro dT 



or 



R' = a Ro T' 



( 22 ) 



v = H R' 



(23) 



Finally, to relate recorded voltage fluctuations to 



T ' 




V' 



(24) 



or 



T ' = C V> 

1 c H 



(25) 



where Cjj is the temperature calibration factor in °C/volt. 
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D. SIGNAL PROCESSING 

The data from the four levels was recorded simultaneously, 
in analog form, on a 14 channel tape recorder. 

The signals were continuously monitored with an 
oscilloscope and an 8 channel strip chart recorder to insure 
the quality of the recorded signal and to detect equipment 
failures . 
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IV, DATA PROCESSING AND ANALYSIS 

A, PRELIMINARY ANALYSIS 

Selection of data periods was based on signal evaluations 
of strip chart records, mentioned in section III.D, and on 
evaluations of experimental logs of equipment status and 
ship positioning. Primary consideration was given initially 
to obtaining 18 minute records with temperature and velocity 
signals at all levels. This was later changed to 21 minute 
records due to characteristics of the anolog spectrum analyzer, 

B . PROCEDURES 

Several procedures were adopted to best handle the 
turbulence data. These are outlined in Appendix A in detail. 

All spectral analysis was performed by processing the recorded 
analog signal on a Federal Scientific Model UA-500-1 Spectrum 
Analyzer-Averager shown in Figure 8. 

The averager output was a logarithmic plot of' the power 
spectral density in Volts 2 /Hz, versus frequency f in Hz. The 
plotter was calibrated so that all plots would concide with 
Bone's (1974) digital plots. Thus, a preliminary visual 
comparison was made to verify the spectra prior to computations. 
Figure 9 illustrates a sample of the plot comparison for a 
data period common to Bone's (1974) study and this study. 

Figure 10 shows a typically good graphical plot of spectrum 
versus frequency used in this study. Figure 11 is an example 
of a poor spectral plot. 
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Figure 8. Spectrum analyzer and related equipment. 




LOG in SPECTRUM 




LOG 10 FREQUENCY 



Figure 9. Comparison of digital and spectrum analyzer 
velocity spectra for 1855-1913L, 18 Jan. 74. 
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Figure 10, Typical good analyzer spectral plot. 
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Figure 11, Typical poor analyzer spectral plot, 
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